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ABSTRACT 

We report on intermittent X-ray pulsations with a frequency of 442.36 Hz from the neutron-star 
X-ray binary SAX J1748. 9-2021 in the globular cluster NGC 6440. The pulsations were seen during 
both 2001 and 2005 outbursts of the source, but only intermittently, appearing and disappearing on 
timescales of hundreds of seconds. We find a suggestive relation between the occurrence of type-I 
X-ray bursts and the appearance of the pulsations but the relation is not strict. This behavior is 
very similar to that of the intermittent accreting millisecond X-ray pulsar HETE J1900. 1-2455. The 
reason for the intermittence of the pulsations remains unclear. However it is now evident that a strict 
division between pulsating and non-pulsating does not exist. By studying the Doppler shift of the 
pulsation frequency we determine an orbit with a period of 8.7 hrs and an projected semi major axis of 
0.39 lightsec. The companion star might be a main-sequence or a slightly evolved star with a mass of 
'--^1 Mq. Therefore, SAX J1748. 9-2021 has a longer period and may have a more massive companion 
star than all the other accreting millisecond X-ray pulsars except for Aql X-1. 
Subject headings: binaries: general-pulsars: individual (SAX J1748.9-2021)-stars: neutron 



1. INTRODUCTION 



Ac c reting millisecond pulsars (AMPs, lAlpar et al.l 
Il982t iBacker e t al.lll98l} ) are transient low mass X-ray 
binaries (LMXBs) that show X-ray pulsations during 
their outbursts. A total of nine AMPs out of 100 
non-pulsating LMXBs have been found to date. The 
reason why only this small subgroup of binaries pul- 
sates is still unknown. The first seven AMPs discov- 
ered showed persiste nt X-ray pulsations throughout the 
outbursts. Rcccntlv iKaaret et al.l ()2006l ) discovered the 
AMP HETE J1900. 1-2455, which has remained active 
for more than 2 years^ but showed pulsations only in- 
termittently during t he first ~ 2 months of activity 
(jCallowav et al.|[2007D . From the transient sourc e Aql 
X-1 pulsations were detected (jCasella et al.l l2007l ) only 
for ~ 150 sec out of the ~ 1.3 Msec the source has (so 
far) been observed with the Rossi X-ray Time Explorer 
( RXTE). 

iGavriil et all ()2006l 120071 ) recently reported on the de- 
tection of ^ 442.3 Hz pulsations in an observation of the 
2005 outburst of a transient source in the globular cluster 
(GC) NGC 6440. The pulsations followed a flux decay 
observed at the beginning of the observation and were 
reminis cent of those observ ed during superbursts; how- 
ever, as IGavriil et al.l (|2007( ) suggest, they could also be a 
detecti on from a new inter mittent accreting millisecond 
pulsar. iKaaret et al.l (j2003l ) report a 409.7 Hz burst oscil- 
lation in an X-ray transient (SAX J1748. 9-2021) located 
also in NG C 6440 and this GC ha rbors at least 24 X - 
ray sources (jPoolev et al.l l200^. so' Gavriil et aT (2007) 
concluded that the burst oscillations and the pulsations 
were probably coming from different X-ray transients in 
the same GC. 

^ Email: diego@science.uva.nl ; Astronomical Institute, "An- 
ton Pannekoek", University of Amsterdam, and Center for High 
Energy Astrophysics, Kruislaan 403, 1098 SJ Amsterdam, The 
Netherlands. 

^ At the time of submitting this letter, the source is still active. 



The exact formation mechanisms behind the pulsations 
of these three sources remains unknown. The existence 
of intermittent pulsations with a small duty cycle implies 
that many other apparently non-pulsating LMXBs might 
be pulsating, bridging the gap between the small number 
of AMPs and the large group of non-pulsating LMXBs. 

We are performing a detailed analysis of all RXTE 
archival data of neutron-star LMXBs to se arch for tran- 
sient pulsations in their X-ray flux (see also lCasella et al.l 
|2007( ) . In this Letter we present the results of our search 
on the three X-ray outbursts observed from the globular 
cluster NGC 6440. 

2. THE NEUTRON-STAR TRANSIENT SAX J1748. 9-2021 IN 
NGC 6440 



NG C 6440 is a GC at 8.5 ± 0.4 kpc (|Ortolani et all 
[T99I . Bright X -ray outbursts from a LMXB were 
repor te d in 1971, 1998 2001 a n d 2005 (iMarkert et all 
1975t lin 't Zand et all 119991: IVCTbunt et all I200C' 
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in't Zand et al 



iMarkwardt fc SwankI 



120011: 

2001h from X-ray and optical obser- 
concluded that the 1998 and 2001 outbursts 
the same object, which they designated 



vations 
were from 
SAX J1748.9-2021 



3. OBSERVATIONS, DATA ANALYSIS AND RESULTS 

We used data from the RXTE Proportional 
Co unter Array (PGA, for instrument information 
see iJahoda et al.l I2006D . Up to July, 2007, there were 
27 pointed observations of SAX J1748. 9-2021, each 
covering 1 to 5 consecutive 90-min satellite orbits. 
Usually, an orbit contains between 1 and 5 ksec of 
useful data separated by 1-4 ksec data gaps due to 
Earth occupations and South Atlantic Anomaly pas- 
sages. Adopting a source position (a = 17''48™52''.163, 
5 ^ -20°2l'32".40; J2000 iPoolev et al.ll200l we con- 
verted the 2-60 keV photon arrival times to the Solar 
System barycenter with the FTOOL faxbary, which uses 
the JPL DE-200 ephemeris along with the spacecraft 
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Fig. 1. — Intensity (2.0—16.0 keV) normalized by Crab vs. time of the three outbursts. Gray symbols show the 16-sec averaged intensity 
during the pointed PCA observations. The continuous line shows the ASM light curve. Black marks at the top mark the times of type-I 
X-ray bursts. Black marks at the bottom mark the times when we detect significant pulsations. Years of the outburst is shown in each 
panel. 



TABLE 1 
Timing parameters for. NGC 6440 



Parameter Value 

Orbital period, Po,-6(hr) 8.764(6) hr 

Projected semi major axis, a^sini (lightsec.) 0.39(1) 

Epoeh of 0° mean longitude\ To (MJD/TDB) 52190.047(4) 

Eccentrieity, e < 0.001 

Spin frequency i/q (Hz) 442.361(1) 

Pulsar mass function, (xIO^'^AIq) 2± 4.8 

Minimum companion mass, Mc (Mq) <^0.1 

^The mean longitude at Q° in a circular orbit corresponds 
to the ascending node of the orbit. 

ephemeris and fine clock correcti ons to provide an 
absolute timing accuracy of 5-8 jis (|Rots et alJl2004[ ). 

We performed a Fourier timing analysis using the 
high-time resolution data collected in the Event 
(E_125us_64M_0_ls) and the Good Xenon modes. Power 
spectra were constructed using data segments of 128, 
256 and 512 seconds and with a Nyquist frequency of 
4096 Hz. No background or dead-time corrections were 
made prior to the calculation of the power spectra, but 
all reported rms amplitudes are background corrected; 
deadtime corrections are negligible. 

3.1. Colors, light curves and states 

From the Standard 2 data (jJahoda et al.ll2006l ). we cal- 
culated colors and intensities with a ti me resolution of 16 
seconds and normalized by Crab (e.g 'Altamiran o et all 
[2007). The PCA observations sample three different out- 
bursts (see Fig.[T]). The color-color diagrams show a pat- 
tern (not plotted) typical for atoll sources. The power 
spec tral fits confirm the identification of these states (see 
also'Kaar et et al.ll2003[ ). We looked for kHz QPOs, but 
found none. 

No thermonuclear bursts were d etected in the first 
outburst, sixteen dur ing the second (|Kaaret et al.ll2003l : 
iGallowav et al.|[2"006D and four during the third one. We 
searched for burst oscillations during all bu rsts in the 15- 
4000 H z frequency range but found none. iKaaret et al.l 
([2.003') reported a 4.4cr burst oscillation at ~ 409.7 Hz. 
We find these authors underestimated the number of tri- 
als by a factor of at least 180, as their estimate did not 
take into account the number of X-ray bursts analyzed 
and the fact that a sliding window was used to find the 
maximum power. Moreover, we also found that the dis- 
tribution of powers is not exponential as these authors 
assumed. Taking into account these effects we estimate 
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Fig. 2. — Top: Dynamical power spectrum of observation 60035- 
02-03-00 showing intermittent pulsations (contours). In the light 
curve (line) three X-ray bursts are seen. The pulse frequency drifts 
due to orbital Doppler modulation. The lowest contour plotted 
corresponds to Leahy power 13 and the highest to 55. The contours 
were generated from power spectra for non-overlapping 128 sec 
intervals of data. Bottom: Leahy normalized (Leahy et al. 198j) 
power spectrum of 512 sec of data centered ~ 7 kscc after the start 
of this observation. Maximum Leahy power is 102, corresponding 
to a single-trial probability of ~ 7 10~ ^^ given P oissonian statistics 
in the photon arrival times ( van dcr Klg[T995l ). Inset: The 2-60 
keV light curve folded at the 2.26-ms period. Two cycles are plotted 
for clarity. The pulse profile is sinusoidal, with a 95% upper limit 
of 0.4% (rms) on the amplitude of the second harmonic. 

the significance for their detection to be < 2.5(7. 

3.2. Pulsations 

We inspected each power spectrum for significant fea- 
tures. We found several, at frequencies ~ 442.3 Hz 
in 7 observations: 60035-02-02-04/05/06, 60035-02-03- 
00/02/03 during the second outbu rst and 91050-03-()7 -00 
during the third outburst (see also lGavriil et al.ll2007l for 
a deta iled analysis of this observation^, lin't Zand et al.l 
(|200lD concluded that the 1998 and 2001 outbursts from 
the LMXB in NGC 6440 were from the same source (Sec- 
tion 2). Since pulsations are detected in both the 2001 
and 2005 outbursts, we can now conclude that these two 
outbursts are also from the same source. Hence, all out- 
bursts observed from NGC 6440 over the last decade are 
from SAX J1748.9-2021. 

The pulsations are detected intermittently, appearing 
and disappearing on time scales of hundreds of seconds. 
The appearance of pulsations seems to be related to the 
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Fig. 3. — Pulse frequency as a function of orbital phase. The 
plot has been obtained by folding all data between the first and 
last pulse detection in 2001. Pulsations were detected during 6 of 
the 18 orbital cycles covered. Drawn curve is the best-fit orbital 
model, measured frequencies and post-fit residuals are shown. The 
residuals' r.m.s. is 1.2 X 10"^ Hz. 

occurrence of type-I X-ray bursts, but the relation is not 
strict. The first two bursts were observed in an observa- 
tion on October 8*'' 2001; the first pulsations a day later. 
During the third outburst we detect four bursts; pulsa- 
tions were only detected after the third one. We also de- 
tected pulsations with no preceding burst. The structure 
of our data does not allow us to tell if pulsations and/or 
other bursts occurred during data gaps. Figure [2] (top) 
illustrates the relation between pulsations and bursts. 
The amplitude of the pulsations varies strongly between 
~ 2% and (often) undetectable (0.3% rms amplitude up- 
per limit at the 95% confidence level). Pulsations are 
seen right after the occurrence of the first and the third 
burst, but in the middle of the data pulsations are present 
without the detection of a preceding burst (although a 
burst could have happened just before the start of this 
data segment). In Figure [2] (bottom) we show a power 
spectrum and corresponding 2-60 keV pulse profile (in- 
set). In these data the pulsation is relatively hard; the 
rms amplitude increases with energy from ^ 1% at 3 keV 
to ~ 3% at 13 keV. 

The 2-10 keV luminosity during the observations in 
which we detected pulsations was between 3 and 4 x 
Iq37 gj.^g g-i (assum ing a fixed Nh = 8.2 x 10^^ cm~^; 
lin't Zand e t al. 2001). Other observations at similar flux 
and those at higher ( up to ~ 5 x 10^^ ergs s~^ in ob- 
servation 91050-03-06-00) and lower fluxes do not show 
pulsations (see Fig.[T]). From 16, 32, 64 and 128 sec aver- 
age colors we found no significant changes in the energy 
spectra correlated with the pulse-amplitude variations. 

We studied the pulse frequency drifts using power spec- 
tra of 128, 256 and 512 sec data and find a clear 8.7 
hours sinusoidal modulation which we interpret as due 
to Doppler shifts by binary orbital motion with that pe- 
riod. In order to obtain an orbital solution, we performed 
a scan on the orbital pararn eters using the method de - 
scribed bv lKirsch et"al1 (|2004D and lPapitto et all (|2005D . 
Our best estimates are listed in Table [TJ The combi- 
nation of data gaps and intermittency of the pulsations 
yielded aliases, which are taken into account by the re- 
ported errors. In Figure [3] we plot the pulse frequency as 
a function of orbital phase. 
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Fig. 4. — Mass-radius relationship for a Roche lobe-filling 
companion (continuous line), isochrones of 0.01, 8 and 12 Gyrs 
with solar motallicity (triangles, crosses and squares, respectively 
; Girar di et al.. .2000) and theo retical Zero-age main sequence 
(ZAMS, dashed line: lTout et aLl irgQeV Black circles mark the in- 
clination of the system as estimated by the mass function of this 
system. 

4. DISCUSSION 

We have discovered intermittent pulsations from the 
neutron-star LMXB SAX J1748.9-2021. Pulsations ap- 
pear and disappear on time scales of hundreds of seconds. 
Although we find a suggestive relation between the ap- 
pearance of the pulsations and the occurrence of type-I 
X-ray bursts (the pulsations appearing after a burst), 
the relation is not strict. We find bursts with no subse- 
quent pulsations and pulsations with no preceding burst 
(although a burst could have occurred in the preceding 
data gaps). From the Doppler shifts on the pulsations we 
determine that the system is in a near-circular orbit with 
period of 8.7 hours and projected radius of 0.39 lightsec. 

The stability of the pulsations (after correcting for the 
binary orbit) strongly suggests that the pulsation fre- 
quency reflects the neutron star spin frequency and that 
SAX J1748. 9-2021 is an accreting millisecond X-ray pul- 
sar. The characteristics of the pulsations are reminis- 
cent of the those found in HETE J1900. 1-2455: in both 
sources the pulsations were only intermittently detected 
and a possible relation between burst occurrence and 
pulse amplitude exists. However, there are differences: 
in HETE J1900. 1-2455 the pulsations were only seen 
during the first two months of the outburst and their 
amplitude decreased steadily on timescales of days after 
the bursts whic h mig ht have caused them to reappear 
(|Gallowav et al.|[200l . In SAX J 1748. 9-2021 we find the 
pulsations in the middle of the 2001 and 2005 outbursts 
and not in the beginning. Furthermore, the amplitude 
of the pulsations behaves erratically, switching between 
detection and non-detection on time scales of hundreds 
of seconds. Despite these differences, the behavior of the 
pulsations in both sources is so similar that we consider it 
likely that the same mechanism causes the intermittency 
of the pulsations in both. 

A related system might be Aql X-1 in which a short- 
lived (^150 s) and very rare (duty cycle of 0.03%) episode 
of strong pulsation s at the neutron-sta r spin frequency 
has been detected (jCasella et al.ll200"7l ). In this source, 
no X-rays bursts were seen in the ^ 1400 s before the 
pulsations, making it unlikely that they were triggered 
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by a burst. It is unclear if the pulsations in Aql X-1 
were accretion-driven or due to unusual nuclear burning 
episodes; the same applies to SAX J1748. 9-2021. The 
extreme rarity of the pulsations in Aql X-1 could indi- 
cate that the mechanism behind them is different from 
that responsible for the pulsations in HETE J1900.1- 
2455 and SAX J1748.9-2021. Nevertheless, irrespective 
of the mechanisms behind the pulsations in these three 
sources, it is clear that a strict division between pulsating 
and non-pulsating sources cannot be made anymore. It 
is possible that all sources pulsate occasionally although 
the recurrence times could be very long. 

Assuming a constant dipolar magnetic field, following 
iPsaltis fc Chakrabartvl ( 1999) (i.e., assuming a geometri- 
cally thin disk and neglecting inner disk wind mass loss, 
radiation drag and GR effects) we estimate the magnetic 
field to be i? > 1.3 X 10^ Gauss. This assumes a 10 km 
radius 1.4M0 neutron star and Mmax, the highest ac- 
cretion rate at which pulses are detected, of 0.28 of the 
Eddington critical value as derived from the luminosity 
obse rved at the time using a bolometric flux correction of 
1.4 (|Migliari fc Fendeil l2006l ) . In the standard magnetic 
channeling scenario, the question remains of what causes 
the large variations in pulse amplitude. 

Comparisons between HETE J1900. 1-2455, 
SAX J1748.9-2021 and the other 7 AMSPs can 
provide clues to understand the pulse-strength vari- 
ations. In SAX J1748.9-2021 and Aql X-1 the time 
scales on which the pulse amplitude can fluctuate are 
as short as ~ 10^ s, too sho rt for the prope rties of the 
neutron star core to change (|Gallowav et al.|[2007i ). So, 
these changes must originate in the di sk or the outer 
layers of the neutron star envelope. IGallowav et al.l 
(11007) suggested (for HETE J1900. 1-2455) that the 
accumulation of matter on the su rface burying the 
magnetic field (jCumming et al.ll2001f ) plays a role. Our 
results show that this mechanism probably cannot 
work for SAX J1748. 9-2021, as the pulsations are not 
seen in the beginning of the outbursts, but instead 
~ 3 weeks and ~ 5 weeks after the start of the 2001 
and 2005 outbursts, respectively, so after a considerable 
amount of matter has already accreted. Interestingly, 
we observe pulsations only around a mass accretion 
rate of ~ 2 x 10"^^ Mq/sbc as inferred from the 
X-ray luminosity, not above or below, indicating that 
instantaneous mass accretion rate rather than total 
accreted mass is the important quantity. 

In both HETE J1900. 1-2455 and SAX J1748.9-2021 
the pulsations seem to appear together with bursts al- 
though the exact connection is complex. This suggests 
that surface processes may affect the magnetic field. 



Hydrodynamic flows in the surface layer of the neu- 
tron star may scree n the magnetic field (see review by 
iBhattacharval [2003 , and references within); perhaps vi- 
olent processes like bursts temporarily affect such flows, 
diminishing screening and enhancing the channeling. 

Alternatively, variations in a scattering or screen- 
ing medium may cause the p ulse amplitude modu- 
lation (s ee e.g. discussi o ns in JP salt is &: Chakrabartyl 
(IT999I). iTitarchuk et all (I2002D. iGogjis et al.l 120071. 
Titarc huk et al. 1 120071 iCaselia et al. I l2007r and references 
within). For our results, the properties of such a medium 
should change on timescales of hundreds of seconds; note 
that we did not detect spectral changes associated with 
pulse strength modulation. 

With an orbital period of ^ 8.7 hours, this binary sys- 
tem is clearly not an ultra-compact binary as usually 
found in globular clusters and in fact, SAX J 1748. 9- 
2021 is the AMSP with the longest orbital period 
after A ql X-1, which has an o r bital period of ~ 
19 hrs (iCh evahe r fc Ilovaiskvl [19911 : I Welsh et lill IMifih . 
The mass-radius relation for a low- mass Ro che lobe- 
filling companion in a binary ()Eggletonlll983t ) is Rc = 

0.24 Mill (72/3 (H-q)l/3 ^2/3/(0 (3 ^2/3+ ^0^(1 + ^1/3))^ 

with Phr the orbital period in hours, Afjvs the mass 
of the neutron star, Rc radius of the companion and 
q = Mc/Mpfs, the mass ratio. Given the mass function 
and the orbital period and assuming a I.4M0 neutron 
star, we plot in Figure |5| the mass-radius relationship for 
the companion star. Given that the age of the Globular 
Cluster NGC 6440 is 10 ± 2 Gyrs (Santos fc Piatti 200|) 
and i ts metallicity is approximately solar ( (Ortolani et aTl 
|1994| ). in Figure m we also plot the isochrones for stars 
with ages of 8 and 12 Gyrs and solar metallicity. Stars 
with a Mc < 0.85 M0 cannot fill the Roche lobe while 
stars with Mc > 0.95 would have a radius exceed- 
ing the Roche lobe. This would imply a donor star mass 
of 0.90 ± O.O5M0. However, for masses of 0.95-1.1 M0, 
stars have evolved off the main sequence so binary mass 
transfer can have affected the radius of the donor star, 
which means we cannot firmly exclude masses of 0.95- 
1 . 1 M0 . Therefore a more conservative mass range for 
the donor star is 0.85-1.1 M0. Intriguingly, this requires 
the inclination to be about 9°, which has a < 1% a priori 
probability for an isotropic sample of binary inclinations. 
Of course, this estimate is assuming that SAX J1748.9- 
2021 is in a primordial binary. If a different evolution- 
ary path took place (e.g. dynamical interactions), the 
mass of the co mpanion might be much smaller (see e.g. 
Ivan Zvl et al][2004) . 
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